Enhanced heating of salty water and ice under microwaves: Molecular dynamics study 
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By molecular dynamics simulations, we have studied the enhanced heating process of salty ice 
and water by the electric field of applied microwaves at 2.5 GHz, and those in the range 2.5-10 GHz 
for the frequency dependence. We show that water molecules in salty ice are allowed to rotate in 
response to the microwave electric field to the extent comparable to those in pure water, because 
the molecules in salty ice are loosely tied by hydrogen bonds with adjacent molecules unlike the case 
of rigidly bonded pure ice. The weakening of hydrogen-bonded network of molecules in salty ice is 
mainly caused by the electrostatic effect of salt ions rather than the short-range geometrical (atom 
size) effect of salt ions since the presence of salt ions with small radii results in similar enhanced 
heating. 
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1. Introduction 

Microwaves arc frequently utilized as a highly efficient 
and less carbon-oxide releasing energy source for cooking 
or drying food, sintering metals and ceramics. It is known 
that dielectric materials including pure and salty water 
are heated efficiently by absorbing the electric field en- 
ergy of microwaves [von Hippel 1954; Hobbs et al., 1966; 
Hasted, 1972; Buchner, 1999; Meissner and Wentz, 2004; 
Takei, 2007]. It is also known that metallic powders and 
magnetic materials such as magnetite and titanium ox- 
ides with oxygen defects TiO^-a; (x > 0) are heated by 
the magnetic field of microwaves [Roy et al, 1999; Pee- 
lamedu et al. 2002; Sato et al. 2006; Suzuki et al. 2008; 
Tanaka et al., 2008]. 

In order to examine the mechanism of microwave heat- 
ing of water and ice at the molecular level, we performed 
molecular dynamics simulations [Tanaka and Sato, 2007]. 
There we showed mechanistically that (i) rotational ex- 
citation of electric dipoles of water coupled with irre- 
versible energy relaxation to the translational energy was 
responsible for the microwave heating of water, (ii) pure 
ice was not heated by microwaves because of strong hy- 
drogen bonds between water molecules, (iii) the exper- 
imentally known enhanced heating of salt water by mi- 
crowaves was due to acceleration of salt ions by the mi- 
crowave electric field, which is usually termed as Joule 
heating. We numerically evaluated the time integrated 
Joule heating term of salt ions E • J and the work associ- 
ated with rotation of electric dipoles E • dP / dt by sum- 
ming up the contribution of all atoms or molecules, where 
E is the electric field of microwaves, J = J^i 1i v i i s the 
current carried by Na + and Cl~ ions, and P = ^ s d s is 
the sum of electric dipoles of all water molecules, with qi 
and Vi the charge and velocity of i-th atom, respectively, 
and d s the electric dipole moment of s-th molecule. We 
obtained the relation |E ■ J = 2|E • dP/dt| for dilute 
saline solution of 1 mol% salinity and 10 GHz microwave 
[Tanaka and Sato, 2007]. 



We stated above that pure ice is not heated by mi- 
crowaves. However, we have a daily experience of pro- 
cessing (melting) frozen food in microwave ovens. Our 
motivation of the present work is to examine the heating 
of salty ice as well as salty water by molecular dynamic 
simulation. In this paper we mainly use the microwaves 
with 2.5 GHz, except for 5-10 GHz microwaves for exam- 
ining the frequency dependence of the heating. 

2. Simulation Procedures 

Ice and liquid water are characterized by a completely 
and partially ordered network of H2O molecules, respec- 
tively, mediated by hydrogen bonds connecting adjacent 
two molecules [Matsumoto et al., 2002]. For this reason, 
we adopt an explicit water model with three point charges 
- the SPC/E water model [Toukan and Rahman, 1985; 
Berendsen et al., 1987]. A typical simulation run requires 
the computation time of a few days for roughly 3000 
water molecules on our Opteron-based cluster machine 
(usually four processors are used per job). Although the 
model has reasonable electrostatic properties, for more 
accurate calculation of the energy absorption, the use of 
TIP4P-FQ model with polarization effect [English and 
MacElroy, 2003; English, 2006] may be preferred but with 
small number of water molecules due to large computa- 
tion times. 

We follow the same numerical procedures as those used 
previously [Tanaka and Sato 2007] . We start all the simu- 
lations with the crystal ice and relax it to the desired tem- 
perature before microwaves are applied. This is because 
molecules of initially random orientations take a very 
long time to arrive at the equilibrium with properly or- 
dered hydrogen-bonded network. A water molecule of the 
SPC/E model in Fig.l illustrates that three atoms form a 
rigid triangle H-O-H with partial charges q(= 0.424e) and 
— 2q residing at the hydrogen and oxygen sites, respec- 
tively. The electrostatic forces are calculated for pairs 
of these atoms while the Lennard- Jones forces are calcu- 




FIG. 1: Three-atom rigid-body model for a water molecule 
(SPC/E model) with partial charges on atom sites. 



lated only for the pairs of oxygen atoms. The equations 
of motion are 
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Here, rj and Vj are the position and velocity of i-th atom, 
respectively, m, and are its mass and charge, respec- 
tively, and Tij = \n — Tj| is the distance between i-th 
and j-th atoms. The second term in the bracket on the 
right-hand side of Eq.(2) serves for two atoms not to 
overlap each other, where <7y is the sum of radii of two 
atoms and &y is the Lennard-Jones energy. For water, 
Sww = 0.65 kJ/mol and a ww = 0.317 nm. Salt ions are 
put randomly in space avoiding water molecules and are 
treated as isolated species with the charge qNa — e, the 
diameter 2a n a = 0.26 nm, and the Lennard-Jones energy 
e Na = 0.062 kJ/mol for Na+ and q C i = -e, 2a C i = 0.44 
nm and e C i = 0.446 kJ/mol for Cl~ [AMBER, 2008]; the 
hybrid rule 
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is used when two species i and j are encountering. The 
diameter of CI ion is comparable with the size of the net- 
work cell whose building blocks are six-membered water 
rings. The SHAKE and RATTLE algorithm [Andersen, 
1983] is adopted in time integration of Eq.(l) and (2) to 
keep the rigid triangular shape of the molecule. The time 
step for integration is 1 fs. 

Our simulation system consists of a cubic box of each 
side 4.2nm that contains (14) 3 water molecules when 
salt is not added. The periodic boundary condition is 
imposed under the constant system volume. Thus, the 
Ewald sum for infinite number of periodic image charges 
is taken to correctly account for the electric field in which 
the particle-particle particle-mesh algorithm is adopted 
[Deserno and Holm, 1998]. The real space cutoff distance 
for both the electrostatic and Lennard-Jones forces is 1.0 
nm, and 32 spatial grids are used for the Fourier-space 
calculation of the nonlocal electrostatic forces. 

After an equilibration phase of 100 ps, a spatially ho- 
mogeneous microwave electric field of the form E mw (i) = 





400 




300 


■ 

-- 






200 


CD 


1 


Qf 










0.0 


— 


■'.ii 






_l_ 
1 


0.0 


y h 






—3.0 


•—• 


3.0 






T 


0.0 




-3.0 



~i 1 1 1 1 r 



to 




o.o Hr^^K^t^^^ 



FIG 



2 4 6 8 10 12 14 
time (x 10 2 ps) 
2: The time history of the heating of salty ice, (a) tem- 



perature of the whole system in Kelvin, (b) the dipole heating 
term (solid line), the Joule heating term of all salt (dotted) 
and that of Na ion (dashed), (c) the x-component velocity of 
Na ions, and (d) that of CI ions with vo = 10 m/s, Qo = 10~ 9 
J. The amplitude of the microwave electric field is shown by 
dotted lines in (c) and (d), where the electric field of 2.5 GHz 
is turned on just after the equilibration phase at time of the 
figure. 



Eosinwi x is applied. The wave frequency / = w/27r 
is 2.5 GHz unless otherwise specified. To overcome the 
numerical noises in the simulation, the used wave am- 
plitude is Eq = 2.3 x 10 6 V/cm which is large but its 
associated energy with the electric dipole of a water 
molecule is less than the thermal energy at room tem- 
perature Eod/ksTsoQK — 0.4. This condition assures us 
the safe use of this electric field in our simulations, where 
d = 7.8 x 10 -30 Cm is the dipole moment of the model 
water molecule. 

3. Results and Discussion 

The initial condition of the simulation is prepared by 
generating crystal ice of the I c structure, and putting 
Na and CI ions of 1 mol% concentration (roughly 3 wt%) 
when salty ice or water is prepared. The I c ice is isotropic 
in three directions whose physical properties are quite 
similar to those of the ordinary ice of the 1^ structure 
[Eisenberg and Kauzman, 1969; Matsumoto et a!., 2002]; 
we swap the O-H bonds of water molecules successively 
until an isotropic state has been reached. The initial 
structure is equilibrated for 100 ps at a given tempera- 
ture, and then the microwave electric field of 2.5 GHz 
frequency is applied. 

Figure 2 shows the time history of (a) the tempera- 
ture (starting at 230 K) of the salty ice system, (b) time 
integrated value of the Joule heating and dipole heat- 
ing terms, (c) the average velocity of Na + along the x 
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time (xl0 2 ps) 

FIG. 3: The time history of the heating of salty water, (a) 
temperature of the whole system in Kelvin, (b) the dipole 
heating term (solid line), the Joule heating term of all salt 
(dotted) and that of Na ion (dashed), (c) the x-component 
velocity of Na ions, and (d) that of CI ions with vo = 10 m/s, 
Qo = 1CT 9 J. The amplitude of the microwave electric field is 
shown by dotted lines in (c) and (d), where the electric field 
is turned on just after the equilibration phase at time of the 
figure. 

direction, and (d) that of CI - . It is noted that the tem- 
perature of salty ice increases monotonically after the 
microwave is turned on. The velocities of salt ions in 
panels (c) and (d) remain small compared to those in 
salty water in Fig. 3, until the ice approaches the melt- 
ing temperature. To clarify the heating mechanism, the 
Joule heating term E • J (dotted and dashed lines for the 
total salt and Na + , respectively) and the dipole heating 
term E • dP/dt (solid line) are plotted separately in (b). 
Evidently, the excitation of dipole rotation is the princi- 
pal heating mechanism of salty ice. 

The heating of salty ice is compared with that of salty 
water in Fig. 3. Salty water with 1 mol% salinity is pre- 
pared in the same way as for salty ice, and is melted and 
equilibrated for 100 ps at 300 K. For the salty water, 
the heating rate increases with temperature above 350 
K, as seen in Fig. 3(a) at which the slope of the curve 
becomes steeper. The velocity oscillations of Na and CI 
ions in (c) and (d) also increase in time and are much 
larger than those for salty ice. The large oscillations at 
later times reveal de-trapping of salt ions from the cell of 
water network, which accelerates the heating due to the 
large amplitude of oscillations. A striking difference is 
found in the heating terms in panel (b): the Joule heat- 
ing term much exceeds the dipole heating term, where 
the former is twice larger than the latter as shown pre- 
viously [Tanaka and Sato, 2007]. It is interesting that CI 
ions with heavier mass and larger radius respond more 
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FIG. 4: The time history of the heating of pure water, (a) 
temperature in Kelvin, (b) the dipole heating term, (c) the x- 
component velocity of hydrogen atoms, and (d) that of oxygen 
atoms with i>i = 1 m/s, Qo = 10 -9 J. The amplitude of the 
microwave electric field is shown by dotted lines in (c) and 
(d), where the electric field of 2.5 GHz is turned on just after 
the equilibration phase at time of the figure. 



than Na ions to the microwave electric field. This im- 
plies that CI ions are not stably contained in the cells of 
the hydrogen-bonded network of water molecules, which 
makes them more easily be de-trapped and move through 
the network. 

As for reference, the heating of pure water against 2.5 
GHz microwave is shown in Fig. 4. In this case, the 
excitation of dipole rotation is the only mechanism of 
energy absorption from the microwave. The heating rate 
of salty ice with 1 mol% salinity is comparable to that 
for pure water. The velocity oscillations are shown for (c) 
hydrogen atoms and (d) oxygen atoms, which are small 
compared with those in Fig. 3 and arc out of phase with 
each other since the center of mass of water molecules is 
fixed because of the charge neutrality. 

As pure ice is not heated by microwaves of the GHz 
band, the presence of salt in ice should be playing a de- 
cisive role in the heating (melting) process in the mi- 
crowave electric field. We showed on the energy basis in 
Fig. 2 that not the Joule heating but the excitation of 
dipole rotation is responsible for the heating. Then it fol- 
lows that electric dipoles of water molecules are allowed 
to rotate in salty ice, which occurs because of the local 
breaking of water network by salt ions. This is well rec- 
ognized in Fig. 5 which shows the time history of the sum 
of the electric dipoles for (a) pure water and (b) salty ice. 
For pure water, molecules begin rotational motions im- 
mediately with the microwave electric field. In the early 
phase when microwaves are applied to salty ice, the ro- 




time (xl0 2 ps) 

FIG. 5: Time history of the sum of all electric dipoles in the x 
direction for (a) pure water (corresponding to Fig.4) and (b) 
salty ice (Fig. 2). Solid lines show the value of electric dipoles, 
and dashed lines show the amplitude of the electric field of 
microwaves. 



(a) (b) 




FIG. 6: 3D view of water molecules and salt ions for salty 
ice (230 K) at (a) the initial stage, and (b) after 100 ps with- 
out microwave application. Chained gray and white spheres 
represent oxygen and hydrogen atoms, respectively, Na and 
CI ions are shown by isolated small and large sphere, respec- 
tively. 



tational motion of the dipoles begins which is somewhat 
retarded and small compared to the case of pure water. 
This is completely different from the case of pure ice. 

Figure 6 shows that the rigid network of ice is distorted 
quickly in the presence of salt ions even before microwave 
is applied, as seen in (b). Experimentally, anomalously 
high permittivity was reported for salty ice [Zhong et al, 
1988]. Here, it is confirmed by numerical simulation at 
the molecular level that the high response of salty ice 
to microwaves is caused by the presence of salt ions in 
ice which weaken the hydrogen-bonded rigid network of 
water molecules. 

The weakening of the hydrogen-bonded network of wa- 
ter molecules which occurs in salty ice is also confirmed 
by specially designed simulation: the diameter of charged 
salt ions is artificially reduced to 0.2 nm so that they 
nicely fit in the cell of the network. Even in such a con- 
dition, heating of salty ice does take place whose heating 
rate is almost three quarters that of the normal case. 
This proves that the weakening of the network in salty 
ice is caused mainly by the electrostatic effect of salt 
ions which is long-ranged, rather than the short-range 
geometrical effect. 
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FIG. 7: The dependence of heating rate on the microwave 
frequency. Open circles and triangels correspond to salty wa- 
ter and salty ice with 1 mol% salinity, respectively, and filled 
circles correspond to pure water at 300 K. 

The dependence of the heating rate on the frequency of 
microwave is examined for pure and salty water and salty 
ice, as depicted in Fig. 7. Microwaves of three frequen- 
cies 2.5, 5 and 10 GHz are used, and the temperature 
increase in the 1.4 ns interval is plotted for pure water 
(filled circles), salty water (open circles) and salty ice 
(triangles). The salinity of salty ice and water is 1 mol%, 
and initial temperature is 300 K for water and 230 K for 
ice. The line fitting yields the frequency dependence of 
dT / dt oc uj 1a and w 1 ' 2 for pure water and salty ice, re- 
spectively. For salty water, we obtain dT/dt oc uj°- 7 . The 
proximity of the scaling laws for pure water and salty 
ice indicates that water molecules in salty ice under mi- 
crowaves rotate as freely as those in pure water. 

4. Summary 

By molecular dynamics simulation, we studied the 
heating process of salty ice, which we daily encounter 
when we melt frozen food in a microwave oven. We 
confirmed that the heating of salty ice is attributed to 
the weakening of the hydrogen-bonded network of water 
molecules in the presence of charged salt ions. In this 
process, the electrostatic effect was dominating over the 
short-range geometrical effect. The dependence of the 
heating rate on the microwave frequency was examined 
for salty ice and water, which again indicated that the 
internal structure of water molecules in salty ice is de- 
formed. Thus, the water molecules in salty ice can rotate 
as freely as those in pure water under the microwave elec- 
tric field. 
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